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^ ■ ABSTRACT 

Q Following the detection of a gravitational microlensing high magnification event 
(HME) in Q2237+0305A attempts have been made to place limits on the dimen- 
sions of the quasar continuum source. The analyses have studied either the observed 
event magnitude or the event duration. The latter approach has been hampered by 
lack of knowledge about the transverse velocity of the lensing galaxy. We obtain both 
upper and lower statistical limits on the size of the continuum source from the ob- 
served HME using determinations of transverse velocity obtained from the published 
monitoring data. Our calculations take account of the caustic orientation as well as 
the component of the caustic velocity that results from stellar proper motions. Our 
determination of source size relies on an estimated duration of 52 days for the HME, 
and so will be refined when more HMEs are observed. We find that the upper and 
lower limits on the magnified region of the R-band continuum source are 6 x 10^^ 
and 2 x 10^^ cm respectively (99% confidence). Through consideration of the joint 
probability for source size and mean microlens mass we find that the mean mass lies 
between - O.OIMq and IMg (95% confidence). 
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j_j ■ The object Q2237-I-0305 comprises a source quasar at red- 
' shift z = 1.695 that is gravitationally lensed by a foreground 
galaxy {z = 0.0394) producing 4 images with separations of 
^1". Each of the 4 images are observed through the bulge 
of the galaxy, which has an optical depth in stars that is of 
order unity (eg. Kent & Falco 1988; Schneider et al. 1988; 
Schmidt, Webster & Lewis 1998). In addition, the proximity 
of the lensing galaxy means that the projected transverse ve- 
locity at the source is an order of magnitude higher than for 
a typical lens configuration. The combination of these con- 
siderations make Q2237-I-0305 the ideal object from which 
to study microlensing. Indeed, Q2237-I-0305 is the only ob- 
ject in which cosmological microlensing has been confirmed 
(Irwin et al. 1989; Corrigan et al. 1991). In 1988 a high mag- 
nification event (HME) was observed in image A. It had a 
measured rise time of ~ 26 days, and was followed by an 
event having a decline time of ^ 3 months. The resulting 
light-curve has been interpreted as a double peaked event 
corresponding to the source having crossed two fold caustics 
inside a cusp (Witt & Mao 1994), implying that the peak 
was caused by a single caustic crossing. The observation thus 
provides an estimate of the event crossing time, and a lower 



limit to the event amplitude. Attempts have been made on 
both of these fronts to interpret the HME in terms of the 
size of the magnified continuum region (eg. Wambsganss, 
Paczynski & Schneider (1990); Ranch & Blandford 1991; 
Jaroszynski, Wambsganss & Paczynski 1992; Webster et al. 
1991). 

The amplitude of an HME is a function of source size 
- smaller sources produce events of larger amplitude. How- 
ever, the event amplitude is also dependent on the width 
or strength of a caustic, a quantity that is the constant of 
proportionality in the near caustic approximation of Chang 
& Refsdal (1979). Witt (1990) termed this constant the flux 
factor and found that it takes on a range of values in the 
complex caustic structures that are produced by high opti- 
cal depth microlensing models. The mean value of the flux 
factor depends on the mass function. Thus the event ampli- 
tudes have a distribution that is model dependent. Wambs- 
ganss, Paczynski & Schneider (1990) found that a microlens- 
ing model with a mean microlens mass of O.225M0 repro- 
duces the observed luminosity variation when the source has 
a radius smaller than ~ 2 x lO^^cm. Several attempts have 
been made to discuss the observed HME amplitude in terms 
of physical models for the continuum region (eg. Ranch 
& Blandford 1991; Jaroszynski, Wambsganss & Paczynski 
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1992). The results of these calculations remain controversial. 
In addition, Lewis et al. (1998) employ the observed differ- 
ence in spectra obtained at two different epochs to infer a 
continuum region that is ^ 10^^ cm (based on a blackbody 
accretion disk model). 

A different approach is to use the observed event du- 
ration in combination with the transverse velocity to ob- 
tain an estimate of source size. Approximate values have 
been obtained assuming a galactic transverse velocity of 
vt ~ 600 fcm see" ^ (providing the caustic velocity with re- 
spect to the source), and a perpendicular caustic crossing 
(eg. Irwin et al 1989; Webster et al. 1991). These analyses 
have yielded values for the continuum source size that are 
consistent with those obtained through consideration of the 
event amplitude. 

In this work we present an analysis of the source size re- 
quired to re-produce the observed event duration from model 
light-curves. This work follows results for the probability 
of the galactic transverse velocity that we have obtained 
previously (Wyithe, Webster & Turner 1999b,c, hereafter 
WWTb, WWTc). In this analysis we consider the effect on 
the caustic crossing time of a stellar velocity dispersion and 
the orientation of caustics with respect to the source trajec- 
tory. 



2 THE MICROLENSING MODEL 

This work is based on an extensive suite of numerical mi- 
crolensing simulations presented previously (WWTc); their 
salient features are discussed below. We adopt the optical 
depth and shear (k and 7) parameters for the lensing galaxy 
of Q2237+0305 calculated by Schmidt, Webster & Lewis 
(1998), and limit our attention to microlensing models in 
which all the point masses have identical mass. Four classes 
of models are produced, combining 0% and 50% smoothly 
distributed matter with source trajectory directions parallel 
to the A-B (7A,7s < 0) and C-D (7a,7s > 0) axes. We 
simulate mock data sets using the sampling rate and period 
of the available monitoring data (Irwin et al. 1989; Corri- 
gan et al. 1991; 0stensen et al. 1995). Two assumptions are 
made for the photometric error: (Ts_b=0.01 mag ((7l_b=0.02 
mag) in images A/B and (7s_b=0.02 mag ((t_l_e=0.04) in im- 
ages C/D. We assume that microlensing is produced through 
the combination of a galactic transverse velocity with an 
isotropic Gaussian velocity dispersion (a* ~ 165 fcm see" ^ 
from Schmidt, Webster & Lewis 1998). Where required, a 
cosmology having S7 = 1 with Hq — 75 km sec~^ is assumed. 

This paper builds on work presented in WWTc. Dis- 
cussions of the microlensing models, the effective transverse 
velocity (Ve//), and the method for the determination of its 
probability Pv{Vgai < Veff\{m)) where (m) is the mean mi- 
crolens mass can be found in that paper. The galactic trans- 
verse velocity is obtained from the effective transverse ve- 
locity by comparing the microlensing rates using histograms 
of light curve derivatives (see WWTb) . 



3 SOURCE SIZE 




Figure 1. A schematic of the geometry of a caustic crossing a 
source. 



3.1 Source Size determination 

The HME in image A observed during late 1988 is described 
by 5 observations (Irwin et al. 1989; Corrigan et al. 1991) 
that constrain a peak in the R-band light-curve. We assume 
that the length of the HME is twice the estimated rise time 
of 26 days (Corrigan et al. 1991). The likelihood for the 
true R-band continuum source size is determined by com- 
paring the observed event duration to distributions of model 
event durations. The calculations of these distributions in- 
clude consideration of the probability for effective transverse 
velocity (WWTc), the normal caustic velocities that arise 
from the proper motion of stars, and caustic orientations. 

The velocity of the caustic network resulting from the 
stream motion of a starfield was calculated by Schramm et 
al. (1993). This result was generalised by Kundic, Witt & 
Chang (1993) to compute the velocity of individual caustic 
points due to the proper motion. We use the expressions 
obtained to find the velocities of caustic points in the yi and 
y2 directions. The vector addition of these components is in 
general in a direction other than normal to the caustic at 
that point. Therefore we also calculate the caustic gradient. 
The method of calculating both the caustic location and 
gradient are described below. 

We implement the contouring algorithm described by 
Lewis et al. (1993). The image curves are followed in discrete 
steps in one direction. When an image line crosses a critical 
curve, the image magnification changes sign due to the op- 
posite parity of critical images on either side of the critical 
curve. The location of the critical curve is then found by 
searching for the image position that produces a lens map- 
ping Jacobian determinant of zero, and used to calculate 
the caustic velocity. The gradient of the caustic tangent is 
found by first computing the location of an image position 
a short distance along the critical curve in both directions. 
The gradient of the caustic tangent is computed from the 
corresponding caustic points. 

In each of the models considered we found the dis- 
tributions of caustic velocities and gradients by searching 
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Figure 2. The probability distribution functions of normal caustic velocities. The thin, medium and thick lines correspond to transverse 
velocities of vt = 0,2.5 and 5x the average projected speed of the stellar velocity dispersion respectively. The thin light line shows the 
Gaussian distribution of stellar speeds. Left: The applied shear is positive (7 > 0). Right: The applied shear is negative (7 < 0). 



for caustics along an equal length of source track in two 
orthogonal directions. The source tracks were each IGrjo 
(corresponding to ^ 70 years for a transverse velocity of 
600 km sec~^) in length, where r/o is the Einstein radius of a 
IMq star. We computed the properties of the caustics that 
intersect with 500 such source tracks in directions parallel 
(7 > 0) and perpendicular (7 < ) to the shear. A separate 
model starfield was used for each track. This is equivalent 
to finding the caustic locations on a 500 x 500 line grid cov- 
ering a IOt^o X 10»7o square of source plane, although the 
use of 1000 starfields makes the sample more statistically 
representative. 

Fig. ^ is a schematic representation of a caustic cross- 
ing, all velocities in this figure are source plane velocities. 
The caustic velocity Vcaust is the vector addition of the 
contribution to caustic velocity that results from stellar 
proper motion (wprop 1, i'piop2), and the transverse velocity 



{Vtranl,Vtran2), where 
(^tran 1 ; '^t?'ari 2 ) — "^trat 



(1) 



Vtran is the galactic transverse velocity relative to the ob- 
server source line of sight, and Ds and Do are the Angular 
diameter distances of the source and lens respectively. The 
normal caustic velocity v±_ is computed from Vcaust and the 
relative orientation between the tangent to the caustic and 
the caustic velocity vector. The magnitude of the normal 
caustic velocity is: 



sin(<^). 



(2) 



where : < </> < -I is the angle between the caustic tan- 
gent and velocity vector. A caustic with a velocity vector 
that makes a small angle with the caustic tangent sweeps 
out a smaller area per unit length per unit velocity than one 
having a velocity normal to its tangent. Therefore, calcu- 
lations of the distribution of normal caustic velocities at a 
stationary source position, include an additional weighting 
of sin(0). 

Figure |^ shows probability distributions of normal caus- 
tic velocities Pc{v±\{m) ,vtran) where vtran — \vtran\- The 
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Figure 3. Bottom: A sample light-curve for image A (7^ = 
-1-0.4) of Q2237+0305. One half of a double peaked event is high- 
lighted by a box. Centre: A top-hat intensity profile (dark line), 
and a Gaussian intensity profile (light line). Top: A plot of the 
the highlighted HME showing the details of the event shapes that 
result from the two source intensity profiles. 
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Table 1. Table showing the upper and lower 99%, 95%, and 90% limits as well as the mean of the distributions ps{Ss |0.1AfQ) and 
Ps{Ss I l.OM©). The values are given for the two directions assumed for the transverse velocity as well as for simulations that contain 
0% and 50% smoothly distributed matter. The limits shown correspond to the cases where the error was crsE = 0.01 mags in images 
A/B, asE = 0.02 mags in images C/D, and where the error was cti^e = 0.02 mags in images A/B, aj^E = 0.04 mags in images C/D (in 
parentheses) . 



0% Smooth Matter 


Mean 
Mass 


Trajectory 
Direction 


S'i„™(99%) 
(xlOl^cm) 


S!o„(95%) 
(xlQl^cm) 


(xlO^^cm) 


Mean 
(xl0i-*cm) 


5„p(90%) 
(x 10" cm) 


5„p(95%) 
(xl0i*cm) 


Sup{99%) 
(xlQi-* cm) 


O.1M0 
I.OMq 
O.1M0 
1.OM0 


7A,7s > 
7A,7s > 
7C.7D > 
7C>7d > 


0.244 (0.238) 
0.394 (0.280 ) 
0.240 (0.237) 
0.370 (0.272) 


0.637 (0.617) 
1.10 (0.746) 

0.622 (0.614) 
1.03 (0.723) 


1.01 (0.972) 
1.80 (1.19) 
0.982 (0.968) 
1.67 (1.16) 


6.70 (6.52) 
12.7 (8.48) 
6.56 (6.50) 
11.9 (8.11) 


12.1 (11.9) 
24.1 (16.0) 
11.9 (11.9) 
22.6 (15.2) 


16.6 (16.2) 
31.3 (21.7) 
16.3 (16.2) 
29.2 (20.5) 


29.2 (29.0) 
49.2 (36.6) 
29.1 (29.0) 
46.0 (35.1) 



50% Smooth Matter 



Mean 
Mass 


Trajectory 
Direction 


Slo^u{99%) 

(xlQl-^cm) 


S';ou,(95%) 
(xlQl-^cm) 


S';o^„(90%) 
(xlO^cm) 


Mean 
(xlO^cm) 


5„p(90%) 
(x 10" cm) 


5„p(95%) 
(xl0"cm) 


5„p(99%) 
(x 10" cm) 


O.1M0 
1.OM0 
O.1M0 
1.OM0 


7a,7s > 
7A,7b > 
IClD > 
IClD > 


0.207 (0.187) 
0.487 (0.274) 
0.206 (0.188) 
0.468 (0.272) 


0.556 (0.493) 
1.45 (0.772) 

0.553 (0.495) 
1.36 (0.760) 


0.895 (0.795) 
2.48 (1.26) 

0.888 (0.796) 
2.26 (1.23) 


5.95 (5.31) 

16.4 (10.1) 
5.73 (5.27) 

14.5 (9.25) 


11.1 (9.75) 
30.1 (19.4) 
10.5 (9.57) 
26.8 (17.7) 


14.6 (13.3) 
38.1 (25.7) 
13.1 (13.2) 

33.7 (23.2) 


24.0 (23.2) 
58.4 (40.4) 
23.3 (23.0) 
49.3 3(6.7) 



left and right hand diagrams correspond to the cases of pos- 
itive and negative galactic shear respectively. Distributions 
are shown in each case for the normal caustic velocities due 
to the stellar proper motions alone (thick line) , and in com- 
bination with a transverse velocity in the yi direction of 2.5 
(medium line) and 5.0 (thin line). Also plotted in Fig. |^ is the 
distribution of stellar speeds (hght line). All velocities are in 
units of the projected mean speed of the isotropic Gaussian 
stellar velocity dispersion. These distributions demonstrate 
the dependence of the normal caustic velocity on the direc- 
tion of transverse velocity; a transverse velocity along the 
direction of caustic clustering (7 < 0) results in lower nor- 
mal caustic velocities because the tangent to the caustic is 
more likely to lie along the source trajectory. Note however 
that the distribution must be independent of the shear di- 
rection if the transverse velocity is zero. 

We assume a circular projected source profile so that the 
event timescale is not dependent on the source orientation. 
We define the source diameter Sa- For a small source {Ss <^ 
rjo) the crossing time of the caustic, and therefore the event 
timescale is 

ATr^Ss/v^. (3) 

Thus, Sa can be estimated by combining a measured 
timescale with statistical information about caustic orien- 
tations and velocities. Eqn ^ demonstrates a systematic un- 
certainty in our determination of source size that is pro- 
portional to the uncertainty in the estimate of the event 
duration. 

Eqn 1^ defines our measure of source size, however its 
interpretation is dependent on the intensity profile of the 
source. The lower panel in Fig. |^ shows an example of an 
extended source light-curve for image A (7A = -1-0.4) com- 
puted using the method described in Wyithe & Webster 
(1999). One half of a double peaked event is highlighted by 
a box and reproduced in the top panel of Fig. H for two 
source intensity profiles: a Gaussian profile (hght line) and 



a top-hat profile (dark line). Cross-sections of the intensity 
profiles are shown in the central panel of Fig. ^ and com- 
parison of these with the corresponding HMEs demonstrates 
that Sa is approximately the source diameter for a top-hat 
source profile, but twice the half width Ss for a Gaussian 
profile. 

Probability distributions for normal caustic velocity 
Pc{v± I (ni) ,Vtran) sxe computcd assuming various values for 
the galactic transverse velocity. By differentiating the cu- 
mulative distribution Pv{V < VtranKm)) (following a cor- 
rection for the contribution of stellar proper motion to the 
value of effective transverse velocity (WWTb)), we also find 
the probability distribution for galactic transverse velocity, 

Pv{vtran = ^^Vtran | (m)) = ^^Pv {Vtran\{m)) . (4) 

The probability distribution for source size Sa is 

Pa{Sa = vj_AT\{m)) = 

tran ■ (5) 

3.2 Source size limits 

Our results are presented in Tab. Values are given for the 
upper and lower 99%, 95%, and 90% limits as well as the 
mean of the distributions Pa{Ss \ O.IMq) and Pa{Sa \ I.OMq) 
for models having various combinations of assumptions for 
the transverse velocity directions and smooth matter frac- 
tion. Limits are shown assuming the small (SE) and large 
(LE, in parentheses) errors defined in Sec. 0. The differential 
distributions of source size Pa{Sa lO.lM©) (thin lines) and 
Ps{Sa I I.OA/0) (thick lines) are plotted in Fig. § The solid 
and dotted lines correspond to the cases of small and large 
errors respectively. We find that the source size lies between 
the values of 2 x lO^^cm and 6 x lO^'^ cm (98% level). 

While the lower source size limit is not very sensitive 
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Figure 4. The differential distributions of source size Ps{Ss \ O.IMq) (thin lines) and ps{Ss \ I.OMq) (thick lines). The solid and dotted 
curves represent the resulting functions when the photometric error was assumed to be asE=0-0^ mag in images A/B and as e=0-02 in 
images C/D, and (J]^e=0.02 mag in images A/B and (7i£;=0.04 in images C/D. 



to the microlens mass assumed, the upper hmit is higher 
for larger assumptions of mean microlens mass. The trend 
is demonstrated in Fig. ^, and is due to the hmit of effec- 
tive transverse velocity being proportional to -y/ (m). Lower 
limits for source size come from the lowest velocity caus- 
tics, which are produced by the stellar proper motions alone, 
hence the insensitivity to the microlens mass. Ifowever the 
upper limits result from the fastest moving caustics, and so 
depend on the microlens mass dependent upper limits to 
transverse velocity. The effect is magnified if a component 
of smooth matter is assumed. Due to the more sparsely dis- 
tributed caustic network produced by such a model, a larger 
transverse velocity is measured (WWTc). Also the contri- 
bution of stellar proper motions is reduced. These factors 
result in an increased sensitivity of source size estimates 
to microlens mass if the model smooth matter content is 
higher. The source size distributions are highly asymmetric. 
Ifowever, due to the stability of the lower end of the distri- 
butions, the mode is consistent, between ~ 1 x 10^'* cm and 
4 X lO^'' cm for all models considered. 

It is apparent from Fig. ^ that the assumptions for 
photometric error and the direction of the transverse veloc- 
ity are important factors, particularly for a larger assumed 



mean microlens mass. Like the assumptions for microlens 
mass, assumptions for the error and transverse velocity di- 
rection effect the measurement of transverse velocity, and 
hence source size. Upper limits are more sensitive than lower 
limits which are set primarily by the stellar velocity disper- 
sion. These effects are also enhanced by the presence of a 
component of smooth matter. 

The co-dependence of mean microlens mass and source 
size is now treated more rigorously. 



3.3 Simultaneous limits of mean microlens mass 
and source size 



In Sec 3.2 the normalised probability distributions, 
Ps(5's|(m)), that the continuum source size is between Ss 
and Ss -I- dSs given that the mean microlens mass is (m) , 
were obtained. The probability that the source size lies in 
the above range and the mean microlens mass is between 
(m) and (m) + d(m) is then given by 



Pm,siSs, (m)) = ps{Ss\{m)) pm{{m)), 



(6) 



where pm{{m}) is the probability that the mean microlens 
mass is between (m) and (m) + d{m). pm{{m}) is a known 
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Figure 5. Contours of percentage peak height of the function pm.s{Ss, {m)). The contours shown are the 3.6%, 14%, 26% and 61% 
levels. The sohd and dotted curves represent the resulting functions when the photometric error was assumed to be asE=0-0^ mag in 
images A/B and crs£;=0.02 in images C/D, and (7^^=0.02 mag in images A/B and o-££;=0.04 in images C/D. The light and dark lines 
correspond to the assumptions of logarithmic and flat priors for transverse velocity referred to in the text. 



function (given in WWTc), that is approximately indepen- 
dent of the source size assumed and was determined under 
the same assumptions for the Bayesian prior (fiat and loga- 
rithmic) of galactic effective transverse velocity as the func- 
tion Pvivtranlim)). However unlike Pv{vtran\{m.)), pm{{m)) 
is dependent on the prior assumed (see WWTc). 

Figure ^ shows contour plots of the two-dimensional dis- 
tribution pm.siSa, (m)) for the assumed source orientations, 
smooth matter fraction and photometric errors. All possi- 
ble combinations of the two values of each parameter are 
included, so that eight models are shown. The contours are 
at 3.6, 14, 26, 61 per cent of the peak height, and so the ex- 
trema of the relevant contours represent 4(t, Sa, 2a, and la 
limits on the single variables. The hght and dark lines corre- 
spond to the assumptions of logarithmic and flat priors for 
effective transverse velocity. The most likely model is that 
(m) ~ 0.05 — O.2M0 (dependent on the model systematics 
of trajectory direction, smooth matter content and photo- 
metric uncertainty) and Sb — 2 x 10^* cm (approximately 
independent of model). 

The most important benefit of combining the distribu- 
tions for Ss and (m) is that an upper limit can be placed 



on the mean microlens mass. The limits on (m) are subject 
to both random and systematic uncertainties; the measured 
mass increases with photometric errors and decreases with 
the fraction of smooth matter. Assuming that the chosen 
values for these parameters bracket the likely models, upper 
and lower limits for (m) can be found by using the most 
extreme values of the contour limits. This yields the 95 per 
cent result that O.OIMq ^ (m) ^ IMq. The upper limit, 
which agrees with that obtained from simultaneous consid- 
eration of probabilities for (m) and Vtran (WWTc) is lower 
than that obtained by Lewis & Irwin (1996), although suffi- 
ciently high that it does not conflict with any popular mod- 
els of galactic halos or bulges. Similarly, the contour maxima 
yield the result that there is a 95 per cent chance of 5 x 10^^ 
cm ^ Ss ^ 10^^ cm. At a given confidence level, both the 
upper and lower limits of source size, and the lower limit 
of mean microlens mass have smaller values in Fig. ^ than 
those obtained from probabilities for single variables in pre- 
vious sections and WWTc. This can be attributed to the 
asymmetry of the distribution in both dimensions. 
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4 CONCLUSION 

The observed peak in the light-curve for image A of 
Q2237-I-0305 in 1988 appears to be one half of a double 
peaked event that is characteristic of a source having passed 
inside of a cusp. If this interpretation is correct then the ob- 
served peak is the result of the source having been crossed by 
a single caustic, and the observations provide a good record 
of the caustic crossing time. The source size is the product 
of the event length and the normal velocity of the caustic 
involved. We have placed limits on normal caustic velocities 
by combining measurements of the galactic transverse veloc- 
ity with calculations of the distributions of caustic gradients 
and the component of the caustic velocity that results from 
stellar proper motions. We find that the most likely size of 
the magnified continuum region is & ~ 2 x 10^"* cm, and 
that the upper and lower limits are 5s 6 x 10^^ cm and 
5*3 ~ 2 X 10^'^ cm (99% level) respectively. Our calculations 
are quite insensitive to the systematic considerations of tra- 
jectory direction, smooth matter content and assumption of 
Bayesian prior for effective transverse velocity (within the 
ranges considered). However, our conclusions are sensitive to 
the microlens mass assumed, as well as to the error assumed 
for the published photometric magnitudes. This sensitivity 
arises from the dependence of the determination of effective 
transverse velocity on these quantities. 

Our conclusions are consistent with those of Wambs- 
ganss, Paczynski & Schneider (1990) and Ranch & Bland- 
ford (1991) who found that the event rise-time and peak 
height require a source smaller than ~ 2 x 10^^ cm (if stel- 
lar mass objects are assumed for the microlenses). The mi- 
crolensing determinations from Q2237-I-0305 support length 
scales of ct ~ 10^* thrs cm associated with X-ray and opti- 
cal changes in continuum emission. Also, theoretical calcu- 
lations of the typical scale-size of a continuum emitting ac- 
cretion disc orbiting a supermassive black-hole (Rees 1984) 
suggest a size of ~ 10^^ cm. 

Simultaneous consideration of microlens mass and 
source size provides an upper limit to the mean microlens 
mass of IA/0 that is unavailable from consideration of the 
microlensing rate alone. We find that O.IM© < (m) < IMq 
(95% level). Unlike the source size, our mass determinations 
are dependent on the prior assumed for effective transverse 
velocity (see Wyithe, Webster & Turner 1999c). 

It was noted in (Wyithe, Webster & Turner 1999b) that 
the statistical determination of the effective transverse ve- 
locity will be aided when monitoring data from a longer 
period becomes available. The narrowing of the probabil- 
ity for the effective transverse velocity will in turn provide 
tighter limits on the mean microlens mass and source size. At 
a known transverse velocity, the component of uncertainty 
due to the range of event lengths that can be produced by a 
source will be reduced upon the observation of more HME 
durations. Thus the extent of the contours in Fig. ^ will be 
reduced as more monitoring data becomes available. 
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